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Abstract 

Recent experiments have observed condensation behavior in a strongly interacting system of 
fermionic atoms. We interpret these observations in terms of a mean- field version of resonance 
superfluidity theory. We find that the objects condensed are not bosonic molecules composed of 
bound fermion pairs, but are rather spatially correlated Cooper pairs whose coherence length is 
comparable to the mean spacing between atoms. We propose experiments that will help to further 
probe these novel pairs. 
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Fermi condensates have been recently observed in dilute atomic gases, first in '^^K h| 

n fi 

and subsequently in ^Li j2|, y]. This new state of ultracold matter represents a Fermi 
gas so strongly interacting that Cooper pairs become correlated in physical space as well 
as in momentum space, similar to the pairs in high-Tc superconductors. Such materials 
are believed to exist in a "crossover" regime, intermediate between weak- coupling (BCS) 
superconductivity and Bose- Einstein condensation (EEC) of tightly bound fermion pairs 0|. 
An ultracold atomic Fermi system is an ideal environment to explore the crossover regime, 
since the effective interactions can be tuned via a magnetic field Feshbach resonance. This 
possibility has led to various predictions based on a "resonance superfluidity" theory of the 
gas 0, ^ El. The BEC limit of the crossover regime was already achieved experimentally 
in the fall of 2003, with the creation of BEC of diatomic molecules composed of fermionic 



atoms 



m 



Because of its close link to high superconductivity, the crossover regime has been a topic 



of intense theoretical investigation, beginning from its prediction 




Il2l| and continuing 
A primary 



through its recent adaptation to ultracold atomic gases 
outcome of crossover theory is that the Cooper pairs begin to become localized in space 
due to many-body correlations as the interparticle interaction becomes large and attractive. 
In the high-Tc superconductor literature, these pairs are referred to as "pre-formed bosons, 
which can exist both above and below their transition temperature to a Bose-condensed 
state. The pairs in the crossover region are smaller than traditionally delocalized Cooper 
pairs, yet are not rigorously bound molecules. In this Letter we explore the link between 
pairs in the crossover regime and molecules by explicitly constructing their wave functions 
for the conditions of the experiment in Ref. We find that the pairs evolve smoothly into 
real molecules as the scattering length is tuned from negative to positive values. We also 
suggest experiments whereby the spatial correlations of the pairs can be probed. Note that 
a recent preprint comes to a similar conclusion for a uniform (i.e., untrapped) Fermi gas 

Q- 

This finding runs counter to the expectations of Refs. llJ], where the pairs are 



identified with actual molecules that are associated with the closed channel wave function 
in two-body scattering theory. If this were the case, then the pair wave function would 
decay exponentially as a function of interparticle separation, regardless of which side of the 
resonance it is on. That this is not the case will be demonstrated below. Additionally, 
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Ref. identified the onset of tlie crossover regime by setting tlie binding energy of the 
molecule, Ti' jma} equal to the Fermi energy Ti'i^'n^nf'l'^ jlm^ where m is the atomic mass 
and n is the number density of atoms. Doing so, one finds that in this regime the scattering 
length (hence the molecular size) is comparable to the interatomic spacing, and the pairs 
are not yet recognizable as distinct molecules. They should rather be considered as spatially 
correlated objects. 

Accordingly, we study in this Letter the correlation length of atom pairs. Our starting 
point is the resonance superfluidity approach [l^ adapted within the Thomas-Fermi de- 
scription 3| • For concreteness, we consider the two- component Fermi gas of ^"i^ near a 
Feshbach resonance between the |9/2 — 9/2 > and |9/2 — 7/2 > states [20|. This system 
possesses a Feshbach resonance whose zero-energy scattering is described by an s-wave scat- 
tering length parametrized by a(-B) = dhgiX " w/ ^B), with a^g = 174ao, w = 7.8 G, and 
AB is the magnetic field detuning in Gauss. For our numerical simulation we have chosen 
the radial frequency z/^ = AOOHz and the trap aspect ratio Ur/^z = 80, as in Ref. Q]. 

The primary objects of the resonance superfluidity theory are the normal and anomalous 
densities, p and k, representing the density of atoms in each species and of correlated pairs, 
respectively, k, can be regarded as the wave function of pairs. Because we work in the local 
density approximation, it is convenient to define these quantities as functions of the location 
R of each pair's center-of-mass in the trap, and the relative momentum k of the pair. The 
equations of motion for these quantities take the BCS-type form 0, tI Ijoj]: 

p(k, R) = n{k, K)u\k, R) + (1 - n(k, K))v\k, R), 



K(k, R) = u(k, K)vik, R)(l - 2n{k, R)) 
u2(k,R)\ 1 ,^ ^ /i(k,R) , 



E(k,R) = VMk,R)2 + A(R)2; ^ , ' =-(l± 



i;2(k,R); 2' E{k,R)' 

h{k, R) = + Vmf{R) + VtrapCR) - Ae (1) 

2m 

/d /r 

where n(k, R) = (exp(£'(k, R)/A;bT) + l)"-*^ is the Fermi-Dirac distribution, Kn/(R) is 
the mean field potential, A(R) is the energy gap, = Airh'^ahg/m, ai,g is the back- 
ground (non- resonant) contact interaction, v = [B — Bo)Afi is the detuning in energy 
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units, g = ^JVhgABAfl is a coupling representing the conversion of free fermions into pairs, 
AB is the field width of the resonance, Afi is the magnetic moment difference between two 
hyper- fine levels of the two-component Fermi gas, Ac is the chemical potential, and Vtrap is 
the external atomic harmonic trapping potential. The chemical potential is fixed by conser- 
vation of the mean number of atoms A^. In this theory, a "molecular field" is introduced 
to simplify the theoretical description of free fermions transforming into spatially correlated 
pairs. Notice that (p is not a distinct physical entity, but is determined once k, is known. 
Thou gh derived for interacting fermions, equations P] can also be applied in the BEC 



limit 
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22j | In this limit it is well- known that (f^(k, R) << 1) and that the BCS 



equation for the gap reduces to the Schrodinger equation for the relative motion of two 
interacting bosons, with energy eigenvalue 2Ac = h'^ /ma^ representing the binding energy of 



the (rigorously bound) bosonic molecules 



21 



After a simple derivation we can extend 



this result to the case of trapped atoms. Then the pair distribution function becomes, in 
the BEC limit, 



mA(R)a2 



iTra 
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where $(R) = y^^A(R) as derived in [2J|. This reference also demonstrated that HFB 
equations transform into the Gross- Pitaevskii equation and that $(R) serves as a solution 
for a molecular BEC interacting through a repulsive interaction with scattering length as = 
2a. This result holds in the BEC limit, where A/|Ac| << 1. (A more careful analysis, 
following 2^, will give corrections, but this is not the main goal of this Letter.) Moreover, 
the first term on the right-hand side of Eqn. (0) is the Fourier transform of the molecular wave 
function —7= — in the relative coordinate r. In this way the same wave function k that 
represents Cooper pairs on the BCS side of the resonance actually represents a condensate 
of real molecular bosons on the BEC side. In general this molecular wave function depends 
on R, meaning that the gas may contain molecules in its high-density center, but correlated 
Fermi pairs at its lower-density periphery. For example in this case twice the chemical 
potential is not quite the molecular binding energy but slightly depends on density Q]- 

The number of pairs can be calculated from the anomalous density as N\, = 
J dUdkuCk, R)^. In the BEC limit K,(k, R)^ ^ p(k, R) for small temperatures, which means 
that almost all atoms are paired. Moreover, using (^El) it is easy to check that in this limit 
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the density of pairs transform into the density of real molecules: 

Pb{R) = J dkK{k,Kf = $(R)2 (3) 

Thus the same function ^(k, R) describes the density of Cooper pairs away from resonance, 
pairs in the crossover regime, and molecules on the BEC side of the resonance as the detuning 
is varied. 

The many- body physics of the crossover regime can be quantified in terms of a 'smooth' 



parameter such as the pair coherence length, usually defined as the rms radius of the pair j21| : 

e'(R) = J drnir, R)VV J drK{r, ^ (A;^(R)/m7rA(R))2 (4) 

Using the above result, it is clear that in the BEC limit < >= in the center of 
the trap (note that the size of a molecule is usually taken instead as the mean value of r. 



< r >= a/2). On the BCS side of the resonance y^^(R) defines the 'size' of the Cooper 
pair. Thus the calculation of the coherence length gives us an insight of how the pairs evolve 
in the crossover regime. 

We present the coherence length versus detuning in Fig. ^ (solid line) for the trap aspect 
ratio, number of atoms, and temperature of the JILA experiment For detunings AB > 
0.5 G the coherence length approaches the familiar BCS result (dash-dot line). For negative 
detunings AB < —1 G on the BEC side of the resonance, the molecular size approaches 
the size evaluated from two-body theory (dashed line). In between, the coherence length 
varies smoothly, illustrating the gradual evolution of pairs into molecules. The sizes of these 
objects remain finite across the resonance, in spite of the divergent behavior of the scattering 
length. This size suppression is the result of many-body physics in the unitarity limit of 
kpa > 1, where the physics of the gas is expected to saturate and to depend only weakly on 
the scattering length. 

To illustrate in greater detail the smoothness of the transition between pairs in the 
crossover regime and molecules, we consider their wave functions, as shown in Fig. |21 This 
figure shows pair wave functions r ■ K,{r, R = 0) in the center of the trap, for detunings 
corresponding to "ordinary" Cooper pairs {AB = LOG, a), pairs in the crossover regime 
AB = 0.1 G, b), and molecules AB = — 0.5G, c). Wave functions of the pairs decay away on 
a length scale set by the coherence length, but in an oscillatory way reminiscent of a damped 
harmonic oscillator. This behavior is a many-body effect, and in fact the scale of the oscil- 
lation is set by the interparticle distance (solid bar). The relative motion of true molecular 



bound states of course decay strictly exponentially, as in Fig. 2c). For small detuning, how- 
ever, the correlation length becomes comparable to the molecular size, and the ringing wave 
functions begin to resemble overdamped oscillators, i.e., they decay exponentially (Fig.2b). 
In this way the character of the pair wave functions evolve smoothly into molecular wave 
functions. This interpretation is somewhat complicated by the fact that the shape of pair 
wave function strongly depends on the trap geometry. Even for small detunings but far from 
the trap's center the coherence length is still large. 

In order to qualitatively understand the JILA experiment Q we now consider the positive 
detuning, BCS side of Fig. [T] We see that at a detuning of around 0.5G the size of the pairs 
becomes comparable to the interparticle distance (dotted line in Fig. Q). This criterion 
marks the crossover regime, where the atom pairs are not momentum-correlated objects 
like BCS Cooper pairs, nor are they yet full-fledged molecules. Significantly, this detuning 
is approximately where a condensate fraction can be observed in the JILA experiment |l|, 
implying that the condensed objects consist of correlated pairs rather then real molecules. To 
estimate the condensate fraction, we assume for simplicity that all pairs are Bose condensed 
at experimental temperatures, so that the condensate fraction is simply Nb/{N/2). The 
true condensate fraction presumably depends on the (unknown) interaction between the 
[airs. This in situ condensate fraction is presented as a function of magnetic field detuning 
in Fig. Q (solid line). This fraction becomes significant only for detunings less than about 
0.5 Gauss from resonance, just where the size of the pairs becomes comparable to the 
interparticle spacing (compare Fig.l). The condensate fraction is quite large near zero 
detuning. In the ideal case of a uniform gas this fraction would be 1 on resonance, but it is 
generally smaller for a trapped gas. 

In the JILA experiment, the Fermi condensate is not directly imaged, but rather is probed 
by a magnetic field sweep that converts the atoms into molecules. This sweep is fast enough 
that it does not affect the many-body properties of the gas, but slow enough that atoms 
are efficiently gathered into molecules. The final detuning is far below resonance (~ — lOG 
in Ref. [l|), so that the molecules are far smaller than the pairs that are being probed. An 
infinitely fast sweep that instantaneously projects pairs onto molecules would therefore not 
yield a significant number of molecules. The condensate fraction, however, could still be a 
significant fraction of unity 17 1. 

In the present calculation, we do not treat the time dependence of the magnetic field, 
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and therefore cannot model the experiment as performed. We can, however, suggest another 
experiment that could probe the crossover regime more fully. Let us consider a hypothetical 
experiment where it would be possible to apply an infinitely fast sweep from a positive 
detuning ABbcs to a final detuning ABbec, i-e., literally projecting pairs onto molecules. 
The final condensation fraction / observed by expansion and imaging will then be defined 
as a product of two probabilities: the first is the projection of the pair wave function ^(r, R) 
onto the molecular wave function (j21), normalized by Nb] the second is the fraction of atoms 
that are paired, {Nb/{N/2)): 

(r 1 \ ^ AT 

JdrdK^ir,R) ^ ^ $(R)j * (5) 

This projection depends not only on the mapping of the wave function of a pair onto the 
wave function of the molecule but also on the condensate wave function as $(R). The 
fermionic condensate wave function /t(r, R) cannot be so easily separated as the product of 
center-of-mass and relative functions as in the BEC case (j2I). It should be said that even 
in BEC case the molecular wave functions in (j2)) will depend on R, which means that the 
molecular size will be different from point to point in the trap. But this dependence is quite 
weak especially for large negative detunings, and is therefore neglected. In the case of a 
large positive detuning ABbcs the size of a Cooper pair is considerably larger than the size 
of a molecule and the number of pairs is quite small itself so the overlap integral (0) will 
be quite small. It is clear that the observed condensate fraction will depend on the geometry 
of the trap as well as on the detuning ABbec of the final point of the sweeping. 

We have calculated the condensate fraction, as seen by this projection technique, for 
three different target molecules defined by ABbec, corresponding to molecules of sizes 
a/2 = 1000ao,500ao and lOOao ( Fig. IH}. As anticipated, the condensate fraction measured 
in this way would be smaller if the pairs are projected onto smaller molecules. Thus there 
are two conditions required to support a large observed condensate fraction: first the ABbcs 
detuning on the BCS side of the resonance should be small enough to support a considerable 
number of pairs compared to the total number of atoms; and second the ABbec on the 
BEC side of the resonance should be chosen so that the corresponding scattering length 
will not be very different from the coherence length corresponding ABbcs ■ The second 
condition means that the size of the pair should be comparable with the size of the molecule. 
For the experiment with "^^K atoms these conditions are fulfilled for ABbcs < 0.6G and 

7 



ABbec > —IG. Of course these results strongly depend on the geometry of the trap and the 
temperature. An experimental map of f {ABbec, ^Bbcs) should prove quite illuminating 
as a probe of the length scales and condensation fractions in the crossover regime. 

In conclusion, we found that the recent experiment 1] can be explained semi- 
quant it at ively by counting the number of Cooper pairs on the BCS side of the Feshbach 
resonance. We suggested a new scenario to probe the crossover regime by mapping the 
condensate of fermionic pairs on the BCS side of the resonance onto molecules on the EEC 
side. We found that a considerable condensate fraction can be observed when the coherence 
length of the pairs is on the order of the interparticle distance. 

This work was supported by the NSF. We acknowledge useful discussions with C. Regal, 
M. Greiner, and D. S. Jin. 
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FIG. 1: Coherence length versus magnetic field detuning (solid line) for K atoms in the JILA 
experiment, in the center of the trap. For comparison, the dashed curve represents the rms molecule 
size a/\/2 corresponding to atoms with a scattering length a = ai,g — The dash-dotted curve 
is the BCS limit of the coherence length. The dotted curve represents the interparticle distance in 
the center of the trap. The trap aspect ratio 1/A = 80, the temperature is T = O.OSTp, and the 
trap contains = 5 x 10^ atoms. 



10 







2x10* 

a.u. 



3x10* 





(C) 



" Ixlo' 2x10" 3x10* 

a.u. 

FIG. 2: Two-body correlation function r • K(r,R = 0) versus interparticle separation r in the 
center of the trap considered in Fig. 1. The panels correspond to the detunings AB = LOG (a), 
AB = O.IG (b), and AB = —0.5 G (c). On the negative detuning side of the resonance, the pairs 
are true molecules. For comparison, the solid bar in (a) shows the interparticle distance in the 
center of the trap. 
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FIG. 3: Condensate fraction versus detuning, as measured by projecting onto molecules, according 
to Eqn. (5), for the trap considered in Fig. 1. The molecules projected onto correspond to atomic 
interactions with scattering lengths 2000, 1000, and 200 (dashed, dot-dashed and dotted lines 
respectively). The solid line is the in situ condensate fraction, given by Nb/{N/2). 
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